Geometric approach to chaos in the classical dynamics of abelian lattice gauge theory 
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A Riemannian geometrization of dynamics is used to study 
chaoticity in the classical Hamiltonian dynamics of a U(l) 
lattice gauge theory. This approach allows one to obtain an- 
alytical estimates of the largest Lyapunov exponent in terms 
of time averages of geometric quantities. These estimates are 
compared with the results of numerical simulations, and turn 
out to be very close to the values extrapolated for very large 
lattice sizes even when the geometric quantities are computed 
using small lattices. The scaling of the Lyapunov exponent A 
with the energy density e is found to be well described by the 
law A tx e 2 . 



I. INTRODUCTION 

Classical dynamical aspects of lattice gauge theories 
have recently attracted some interest @,|. The classical 
limit of a lattice gauge theory is interesting both from 
the point of view of classical dynamical system theory 
and in the perspective of its quantum counterpart. As 
far as the latter aspect is concerned, the interest mainly 
resides in the fact that very few non-perturbative tools 
are available to study quantum gauge field theories, while 
in the classical limit it is in principle possible to exactly 
simulate the real time evolution of the system at any 
energy. From the point of view of the theory of classi- 
cal dynamical systems, lattice gauge theories are highly 
non-trivial many-degree-of-freedom Hamiltonian systems 
which exhibit a rich and interesting phenomenology. The 
classical Hamiltonian dynamics of such systems is known 
to be chaotic [|lj; however, a precise characterization of 
the different chaotic regimes which may occur in these 
systems is still lacking. 

Many particular aspects of this general problem have 
been considered in the literature, concerning the proper- 
ties of pure gauge theories and of theories coupled with 
matter (mainly Higgs fields), e.g., the systematic study 
of the Lyapunov spectra || , the study of thermalization 
processes |J, and the relation between Lyapunov expo- 
nents and observable quantities like the plasmon damping 
rate §]. 



A particular problem which is still open is the depen- 
dence of the largest Lyapunov exponent A of the pure 
Yang-Mills lattice theory on the energy density e — en- 
ergy per plaquette, or energy per degree of freedom — 
particularly at low e HQ. First, the e scaling of A seems 
different according to the fact that the theory is abelian 
- U(i) — or non-abelian — SU(2) and SU(3); in the 
latter case two different scalings have been measured, 
namely A oc e 1 / 4 || and X oc s 0, while in the former 
case a rapid decrease of A at low e was observed 0j . As we 
shall see in the following, our results suggest that in the 
U(l) case the power law A oc e 2 holds. As pointed out by 
Miiller and Trayanov || and subsequently by Nielsen et 
al. 0, such a problem is interesting because it is tightly 
related with the problem of the relevance of the chaotic 
lattice dynamics for the continuum limit of the gauge 
theory. In fact, intrinsic dimensional arguments can be 
used to show that the lattice spacing a, which can be 
set equal to one in all the numerical simulations after a 
convenient choice of units, enters the relation between A 
and e as follows: 



aA(a) = f{ae(a)) , 



(1) 



hence if one observes numerically a power law A oc e k , 
the latter can be read as 



A(a) oc a k 1 e(a) 



(2) 



This means that considering a continuum limit a — > in 
which e(a = 0) is finite, corresponding to a finite temper- 
ature of the resulting field, then the Lyapunov exponent 
is finite in the limit a — > only for the particular expo- 
nent k = 1. Larger exponents k > 1 would imply that 
lim a ^o A (a) = 0, thus leading to a regular continuum 
theory, while exponents k < 1 would mean that in the 
continuum theory the Lyapunov exponent diverges. The 
linear behavior plays then a very special role. As regards 
non-abelian theories, in Ref. Q some evidence is reported 
supporting the fact that the correct scaling is the linear 
one, A oc e, the other one (A oc e 1 / 4 ) being a spurious re- 
sult due to finite-time effects. According to these results 
the continuum limit, for small but finite energy densities, 
of the Lyapunov exponent of non-abelian lattice gauge 
theories is finite. In fact, extracting reliable informations 
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from numerical simulations in the low-energy regime is 
very difficult, mainly because of finite-time effects, which 
become very important at small energy densities as the 
characteristic instability time scales grow, and of finite- 
size effects which are known to be rather large in typical 
lattice gauge theory simulations. 

In the present work we apply a recently proposed for- 
malism |10j-|l4|1, which is based on a Riemannian ge- 
ometrization of Hamiltonian dynamics, to the classical 
dynamics of lattice gauge theories. Such a formalism al- 
lows one to relate chaotic dynamics and curvature prop- 
erties of suitable manifolds, and to obtain an analytic 
formula for the Lyapunov exponent in terms of average 
curvature fluctuations JlJJl^]. The quantities entering 
this formula are statistical averages which can be com- 
puted regardless of the knowledge of the dynamics, either 
by MonteCarlo or molecular dynamics simulations, or an- 
alytically in some cases . As a first step, we apply this 
formalism to the abelian — U(l) — lattice gauge theory, 
which is the simplest one, leaving the non-abelian case to 
future work. In the case of a U(l) gauge theory we per- 
form a precise numerical measurement of the Lyapunov 
exponent by simultaneous integration of the Hamilton's 
equations and of the tangent dynamics equations using a 
precise symplectic algorithm Eqj for several lattice sizes. 
In these simulations we measure also the relevant geo- 
metric observables that allow for a characterization of 
the dynamics and that enter the above-mentioned theo- 
retical expression for the Lyapunov exponent. We find 
that the analytic estimate compares very well with the 
outcomes of the numerical simulations for the Lyapunov 
exponents. Moreover, we find that the theoretical esti- 
mate is almost free from finite-size effects; for small en- 
ergies we find the law A oc s 2 already when inserting in 
the formula the values of the geometric observables ob- 
tained with small (4x4x4) lattices, while the numerical 
values of the Lyapunov exponents are affected by large 
finite-size effects. 

The paper is organized as follows: in Section |l| we 
review very briefly the geometric theory of Hamiltonian 
chaotic dynamics; in Section III we describe the model 
and the observables studied. Section IV is devoted to 



the presentation and the discussion of the results, and in 
Section [v] we draw some conclusions and we outline some 
future developments. 



II. GEOMETRY AND CHAOTIC DYNAMICS 

Let us now recall very briefly the main points about 
the geometric theory of Hamiltonian chaos. Details can 
be found in Refs. |lQ-M. Despite the fact that this the- 
ory is still at its beginning, it has already proved useful 
not only in its original context, but also in connection 
with the problem of the relationship between dynamics 
and statistical mechanics (in particular phase transitions) 

MM- 



Hamiltonian dynamics can be rephrased in geometrical 
terms owing to the fact that the trajectories of a dynam- 
ical system with quadratic kinetic energy can be seen as 
geodesies of a suitable Riemannian manifold. There are 
several choices for the ambient manifold as well as for 
the metric tensor. As already discussed in Ref. |Tl|-[l^] a 
particularly useful ambient space is the enlarged config- 
uration space-time M x IR 2 , i.e. the configuration space 
{q 1 , . . . , q 1 , . . . , q N } with two additional real coordinates 
q° and q N+1 . In the following q° will be identified with 
the time t. For standard Hamiltonians Tt = T + V(q) 
where T = ^aijtfcf , this manifold, equipped with Eiscn- 
hart's metric gE, has a semi- Riemannian (Lorentzian) 
structure (det gE — —1). The arc-length is given by 



ds 2 = aiJ dq l dq J - 2V(q){dq u y + 2dq l) dq 



n . N+i 



(3) 



where both i and j run between 1 and N. Let us restrict 
to geodesies whose arc-length parametrization is affine, 
i.e. ds 2 = c 2 dt 2 ; simple algebra shows that the geodesic 
equations 



1^ 



become Newton equations (without loss of generality 



Sij is considered) 



d 2 q l 
~dt 2 



dV 



(5) 



for i = 1 . . . iV, together with two extra equations for q° 
and q N+1 which can be integrated to yield 



q°=t 

q N+1 = ^t + r, - / I(q.q),// 



(6a) 
(6b) 



where L(q, q) is the Lagrangian, and ci, C2 are real con- 
stants. In the following we set c\ = 1 in order that 
ds 2 = dt 2 on the physical geodesies. As stated by Eisen- 
hart theorem |[L9| , the dynamical trajectories in config- 
uration space are projections on M of the geodesies of 
(M x R 2 ,3 E ). 

In the geometrical framework, the stability of the tra- 
jectories is mapped on the stability of the geodesies, 
hence it can be studied by the Jacobi equation for 
geodesic deviation 



D 2 J 
~d^ 



+ R(ff, J)7 = , 



(7) 



where R is the Riemann curvature tensor, j is the ve- 
locity vector along the reference geodesic j(s), D/ds is 
the covariant derivative and J, which measures the de- 
viation between nearby geodesies, is referred to as the 
Jacobi field. The stability — or instability — of the dy- 
namics, and thus deterministic chaos, originates from the 
curvature properties of the ambient manifold. In local 
coordinates, Eq. ^ is written as 
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D 2 J^ 
ds 2 



„ dq v n dq a 
vp ° ds J ds 



= 



and in the case of Eisenhart metric it simplifies to 
d 2 J l d 2 V 



dt 2 



dq t dq- 



■P =0 , 



(8) 



(9) 



which is nothing but the usual tangent dynamics equa- 
tion for standard Hamiltonians. The Lyapunov expo- 
nents are usually computed evaluating the rate of expo- 
nential growth of J by means of a numerical integration 
of Eq. |j) (20 1 . 

In the particular case of constant curvature manifolds, 
Eq. (Q) becomes very simple 



D 2 J^ 
ds 2 



KJ^ = 



(10) 



and has bounded oscillating solutions J ~ cos(\/K s) or 
exponentially unstable solutions J ~ exp(\/ — K s) ac- 
cording to the sign of the constant sectional curvature 
K, which is given by 



K = 



Kr 



K 



N-l N{N-l) 



(11) 



where K R = R^q^q" is the Ricci curvature and 1Z 
is the scalar curvature; R^ v is the Ricci tensor. Manifolds 
with K < are considered in abstract ergodic theory (see 
e.g. Ref. p2[). Krylov [g3[ originally proposed that the 
presence of some negative curvature could be the mecha- 
nism actually at work to make chaos in physical systems, 
but in realistic cases the curvatures are neither found 
constant nor everywhere negative, and the straightfor- 
ward approach based on Eq. (|h]) does not apply. This 
is the main reason why Krylov's ideas remained confined 
to abstract ergodic theory with few exceptions. 

In spite of these major problems, some approximations 
on Eq. (0) are possible even in the general case. The key 
point is that negative curvatures are not strictly neces- 
sary to make chaos, and that a subtler mechanism related 
to the bumpiness of the ambient manifold is actually at 
work. Upon an assumption of quasi-isotropy of the am- 
bient manifold, i.e., that the manifold can be obtained as 
a small deformation of a constant-curvature space (see 
Ref. |l| for details), Eq. (@) can be approximated by an 
effective scalar equation which reads 



d 2 i\) 
~dt 2 



+ K(t) ip = 



(12) 



where ip is a generic component of the vector J (in this 
approximation all the components are considered equiv- 
alent), and K (i) is a stochastic process which models the 
curvature along the geodesic curve. Such a stochastic 
model is defined by 



K{t) = {k R ) + {8 2 k R ) 1 l 2 r){t) 



(13) 



where k R = K R /N, (•) stands for an average taken along 
a geodesic, which, for systems in thermal equilibrium, 
can be substituted with a statistical average taken with 
respect to a suitable probability measure (e.g. the micro- 
canonical or the canonical measure); rj(t) is a station- 
ary (^-correlated Gaussian stochastic process with zero 
mean and variance equal to one. Using Eisenhart met- 
ric, and for standard Hamiltonians, the non-vanishing 
components of the Riemann tensor are Roioj = 9 qi d qj V, 
hence the Ricci curvature has the remarkably simple form 



k R = 



(14) 



where V 2 is the Euclidean Laplacian operator. Equa- 
tion ( |l2| ) becomes a stochastic differential equation, i.e. 
the evolution equation of a random oscillator p4[ |. It is 
worth noticing that Eq. ([l2|) is no longer dependent on 
the dynamics, since the random process depends only on 
statistical averages. The estimate of the Lyapunov ex- 
ponent A is then obtained through the evolution of the 
second moments of the solution of ( |l2|) as 



A = lim — lo, 

t — >oc 2£ 



(15) 



As shown in Ref. fl 
sion for A: 



^ 2 (0) + V 2 (0) 
, this yields the following expres 



A(fc,er fe ,r) 



1 



A 



4fc 
3A 



where 



A 



<j1t 



64/fc 3 



1/3 



2y/k(k + (Tfe) + 7T(Tfc 



(16) 



(17a) 



(17b) 



in the above expressions k is the average Ricci curvature 
k = (kR) and o~k stands for the mean-square fluctuation 
of the Ricci curvature, at = (S 2 k R ) 1 ^ 2 . 

The advantages in using the geometric approach to 
Hamiltonian chaos are thus evident. In fact, it is possible 
to give reliable estimates of the Lyapunov exponent with- 
out actually computing the time evolution of the system: 
the estimate ( |T6| ) of A depends only on statistical aver- 
ages which can be cither computed analytically in some 
cases (for instance in the case of the FPU model |^2[) or, 
in general, extracted from a Monte Carlo or a dynamical 
simulation, as it is the case of the model studied in the 
present work. 

The behavior of the average geometric observables as 
the control parameter (e.g. the energy density or the 
temperature) is varied conveys an information which goes 
beyond the possibility of computing the Lyapunov expo- 
nent. In fact, one can look at the random oscillator equa- 
tion (jl^) as an effective Jacobi equation for a geodesic 
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flow on a surface M whose Gaussian curvature is given by 
the random process Kit). As long as nonlinear coupled 
oscillators are considered, the average Ricci curvature is 
positive, hence M can be regarded as a sphere with a fluc- 
tuating radius. In the limit of vanishing fluctuations, one 
recovers the bounded evolution of the Jacobi field associ- 
ated with integrable dynamics. Chaos suddenly appears 
as curvature fluctuations are turned on, nevertheless it 
it will be "weak" as long as <Jk <C k, i.e. as long as M 
can be considered as a weakly perturbed sphere. On the 
contrary as the size of curvature fluctuations becomes of 
the same order of the average curvature, o~k ~ k, M can 
no longer resemble a sphere, and the dynamics will no 
longer "feel" the integrable limit. Hence we expect the 
dynamics to be strongly chaotic. This is by no means a 
deep explanation of the existence of weakly and strongly 
chaotic regimes in Hamiltonian dynamics. Nevertheless 
it shows how the simple geometric concepts which en- 
ter the Riemannian description of Hamiltonian chaos, 
besides providing effective computational tools, are also 
useful in helping one's physical intuition with images and 
analogies which would be difficult to find elsewhere. 



We can pass to a standard Hamiltonian formulation by 
putting 



P=^=ag 2 U 
dU 



(21) 



thus obtaining 

H=\Y<( P ^ P ^)+ E (l-^TrlV 

plaqucttcs 



a fj 



links 



(22) 



The parameters a and g 2 can be scaled out, so we set 
g = a = 1. 

The Hamiltonian (p2|) is the classical Hamiltonian for a 
lattice gauge theory with a generic gauge group G. Let us 
now specialize to the abelian group G = U(l). Choosing 
the representation 



U 



cos ip sin ip 
— sin ip co sip 



(23) 



we have 



III. MODEL AND DYNAMICAL OBSERVABLES 

The dynamical system that we are now considering is 
the classical lattice gauge theory based on the abelian 
gauge group U(l). The Hamiltonian of such a system 
can be derived in two ways, either as the dual (Wegner 
model) of a lattice planar spin system pa], or starting 
from the Kogut -Susskind Hamiltonian for a generic gauge 
group G [|25|,[26| and then specializing to the group U(l). 
In the latter case we obtain the same Hamiltonian as 
in the former case by choosing the 50(2) real matrix 
representation of the group U(l). 

The lattice Lagrangian, obtained from the Wilson ac- 
tion by fixing the temporal gauge (all the temporal com- 
ponents of the gauge fields are set to zero) and then by 
taking the continuum limit in the temporal direction, is 



£ = 



ag 

links 



^ plaqucttcs 



(18) 

where U x ^ <E G is a group element defined on a link of 
a (i-dimensional cubic lattice, labeled by the site index x 
and the oriented lattice direction [i, g 2 is a coupling con- 
stant, a is the lattice spacing, (•, •) stands for the scalar 
product between group elements, defined as 



(A,B) = lTr(AB^) 



(19) 



and is a shorthand notation for the plaquette oper- 
ator, 



U^iy — U Xl i±U x -\-ii^i/U- t 



X-\-fl-\-l^, — flU X -\-b'. 



(20) 



P = U = <p[ ~ SmV C ° SV ) (24) 
r 1 — cos ip — sin ip 1 



and we find 



(25) 



To write the plaquette operator, we use the fact that the 
group U(l) ~ 50(2) is abelian. Then the product of the 
four rotations is a rotation of the sum of the oriented 
angles along the plaquette (the discrete curl of the field 

V) 



(26) 



and the magnetic energy, i.e., the potential energy of the 
dynamical system can be written as 



V 



( 1 - ^Tr Up V ) = ^ (1 - cos^) 



plaqucttes 



plaqucttcs 



(27) 



Our canonical variables are then the angles <p x ,ii and the 
angular momenta tt x ^ = (p X4l , and the Hamiltonian J2^) 
becomes 



n = 2 E <m + E C 1 - cos ^.f") ■ ( 28 ) 

links plaqucttes 

Constant-energy (microcanonical) simulations have been 
performed on three-dimensional lattices — with lattice 
sizes ranging from 4 3 to 15 3 — by integrating the canon- 
ical equations of motion 
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fx,, 



'X.fl , 

dV 



(29a) 
(29b) 



where x runs over the L 3 lattice sites and [i = 1,2,3, 
by using a precise third-order bilateral symplectic algo- 
rithm jl5[. We remind that symplectic algorithms are 
integration schemes which exactly preserve the canonical 
form of the equations of motion. The exact Hamilton 
equations are replaced, in the time discretization proce- 
dure, by a map that is symplectic, hence the discrete-time 
flow that approximates the true Hamiltonian flow gener- 
ated by Hamilton' s equations is still Hamiltonian. All 
the geometric constraints on the phase space trajectories 
which are enforced by the canonical form of the equations 
of motion are thus exactly preserved during the numer- 
ical integration procedure and the Hamiltonian flow is 
faithfully represented. As a by-product of this features, 
symplectic algorithms conserve very well the total energy 
of the system: in our simulations relative energy fluctu- 
ations were of the order of 1CP 7 -j- 10~ 8 . In Eqs. ( |29| ) 
V is given by Eq. (p7|), whose explicit expression on a 
three-dimensional lattice is 



(30) 



where (fjtv) = 12, 13, 23. The forces (rhs of Eq. |29bJ) are 
given by 



dV 

dfx,p 



2J Bmtpx-n-s^+s - sxacfx^+s . (31) 



5=1.2 



In order to compute the largest Lyapunov exponent A by 
the standard method |po| , the tangent dynamics equa- 
tions (H), which now reads as 



Jx„ 



££ 

V v 



d 2 V 



Jy.v — 



(32) 



have been integrated simultaneously with the Hamilton's 
equations ( p9| ) by means of the same algorithm. The 
largest Lyapuonov exponent has then been computed ac- 
cording to the definition 



A = lim — log 

t^oo t 



\j\ 2 (t) + \J\ 2 (t) 



1/2 



|jp(o) + |jp( ) 



1/2 



(33) 



Here |J| = J2 X Jx,p. ^ the squared Euclidean norm 
of the tangent vector J. 

According to the discussion of Sec. ||, the relevant 
geometric observable which is able to characterize the 
chaotic dynamics of the model is the Ricci curvature kn, 
computed with the Eisenhart metric, defined by Eq. ( 14 ) , 
which now can be rewritten as 



COB tp XllJa/ + COSifx-u^v 



(34) 



The average and the rms fluctuations of the Ricci curva- 
ture kn have been computed in all the simulations. The 
results are presented and discussed in the next Section. 



IV. RESULTS AND DISCUSSION 

The average and the fluctuations of the Ricci curvature 
(|34| ) are plotted in Fig. [l] against the energy density e 
for three different lattice sizes. The resulting patterns re- 
veal that in the low-energy regime the average k is much 
larger than the fluctuation Uk , whereas in the high-energy 
regime the converse is true. These results are qualita- 
tively the same as those that have been recently obtained 
for a planar (XY) spin model on a 2 — d cubic lattice [jl6| . 
Such a fact is not surprising, since the present model, 
possessing a local gauge invariance under proper planar 
rotations, is expected to behave on a e?-dimensional lat- 
tice — from a statistical-mechanical point of view — as 
its dual model (invariant under a global symmetry) on 
ad— 1-dimensional lattice [ p5| |. Moreover, this means 
that our results are consistent with the fact that the sys- 
tem undergoes a Berezinskij-Kosterlitz-Thouless (BKT) 
transition at a finite energy density (of order e ~ 1 in our 
units). However, the analysis of the geometric quantities 
is not expected to give sharp indications of the presence 
of a transition in the case of a BKT transition. On the 
contrary in the case of second-order transitions the fluc- 
tuations of the Ricci curvature exhibit peculiar cusp-like 
behaviours |R^,[n| which are completely absent here. 

From the arguments reported at the end of Sec. || 
we expect a crossover from weak to strong chaos where 
<7fc ~ k. The results reported in Fig. [j] thus suggest that 
such a crossover should occur around e ~ 1, perhaps at 
a somewhat higher value of e in the case of a 15 3 lat- 
tice. Actually, as it is shown in Fig. |^, in the energy 
region around e ~ 1 the numerically computed values of 
the largest Lyapunov exponent rapidly grow. At larger 
energies the values of A decrease because as e — > oo the 
model becomes integrable. At variance with what is usu- 
ally found in other models (where the X(e) pattern is 
rather stable when the number of degrees of freedom of 
the system is varied) sizeable variations in the A(e) pat- 
tern are observed here at different lattice sizes. Again, 
the phenomenology is very close to the one observed in 
the XY model on a 2 — d lattice. The fact that there is a 
change in the chaotic properties of the dynamics is well 
evident looking at Fig. ^, where the Lyapunov exponents 
are plotted against £ in a log- log scale. In fact, in an en- 
ergy range around e ~ 1, the behaviour of A as a function 
of the energy density deviates from a steep power law to 
a smoother one. In Fig. |^ the numerical values of A 
are compared with the theoretical estimates obtained ac- 
cording to the geometric theory outlined in Sec. [□], i.e., 
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obtained by substituting the computed values of k and 
07c, shown in Fig. Q, into Eq. ©. 

Two facts are immediately evident from Fig. ^. First, 
in the low-energy region (weak chaos) the theoretical esti- 
mates show a power-law behaviour A oc e 2 already using 
the geometric values computed using a very small (4 3 ) 
lattice, and the exponent of the power law remains the 
same for all the lattice sizes, whereas the numerical val- 
ues of A show a less steep e-dependence for small lattices. 
Second, both the low-energy power law behaviour and the 
actual values of the numerical Lyapunov exponents are 
the closer to the theoretical estimates the larger the lat- 
tice size is. Let us mention that the actual value of A is 
theoretically underestimated in the transition region, as 
it is well evident from Fig. [|. Without entering the de- 
tails of this mismatch, we mention that the reason can be 
understood in the light of a very similar situation found 
in 1 - d, 2 - d and 3 - d XY models 

Moreover, the actual values of the theoretical estimates 
are almost free from finite-size effects, which are instead 
very large in the numerical values of the Lyapunov ex- 
ponents. This fact is particulalrly evident at very low 
energy densities, as shown in Fig. where the numer- 
ical and theoretical values of A are plotted against the 
lattice size. 



V. CONCLUDING REMARKS 

The results presented in the previous section clearly 
show that the geometric approach to Hamiltonian chaotic 
dynamics outlined in Sec.n] is particularly efficient to 
provide very good estimates of the largest Lyapunov ex- 
ponent for a classical lattice gauge theory with abelian 
gauge symmetry. In particular, the theoretical estimates 
of the Lyapunov exponents are almost free from finite-size 
effects, especially at low energies. The theoretical values 
are very close to the numerical values that one extrapo- 
lates for very large lattice sizes even when the geometric 
quantities are computed using very small lattices. Also 
the dependence of A on the energy density in the weakly 
chaotic (low-energy) region is A oc e 2 already for a 4 3 
lattice, where the numerical values of the A's exhibit a 
much less steep scaling. 

This is a confirmation of the fact that the geometric 
estimates of the Lyapunov exponents obtained by Eq. 
( |l6| ) are very good for large systems, as already observed 
in other cases p3|. It is not surprising since the whole 
theory uses the large size of the systems as a hypothe- 
Moreover, this feature is very promising in 
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the perspective of applying the geometric theory to non- 
abelian lattice gauge theories. In fact, as already men- 
tioned in the Introduction, in the non-abelian case there 
is still some doubt about the correct value of the expo- 
nent a of the power law A oc e a in the low-energy region, 
and this has remarkable importance for the relevance of 
chaos in the continuum limit of the theory. The applica- 



tion of the geometric theory of chaos might perhaps help 
in obtaining reliable estimates of the exponent a already 
using small lattice sizes. 
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FIG. 1. Plots of the average k (full circles) and of the rms 
fluctuation at (open circles) of the Ricci curvature vs. the 
energy density e for three different lattice sizes: from top to 
bottom, L 3 = 4 3 , 8 3 , 15 3 . Errorbars are smaller than the size 
of the data points. 
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FIG. 2. Numerically computed largest Lyapunov exponent 
A vs. the energy density e for three different lattice sizes: 



FIG. 3. Largest Lyapunov exponent vs. the energy den- 
sity e for three different lattice sizes: from top to bottom, 



from top to bottom, L =4,8' 
than the size of the data points. 



15 3 . Errorbars are smaller L 3 = 4 3 ,8 3 ,15 3 , in a log-log scale. The full circles are the 

numerically measured values, and the open circles are the 
theoretical estimates obtained by inserting the average and 
fluctuations of the Ricci curvature, reported in Fig. [j], in Eq. 
(p"6|). The low-energy behaviour A oc e 2 is well evident from 
the theoretical estimates already for a 4 3 lattice, where the 
scaling with e of the numerically computed Lyapunov expo- 
nents is completely different. 




FIG. 4. Numerical results for the Lyapunov exponent at 
e = 1CP 2 (full circles) compared with the theoretical estimate 
obtained by Eq. ( [lrj ) (open circles) vs. lattice size. 







